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We study binary mixtures of ultra-cold atoms, confined to one dimension in an optical lattice, 
with commensurate densities. Within a Luttinger liquid description, which treats various mixtures 
on equal footing, we derive a system of renormalization group equations at second order, from which 
we determine the rich phase diagrams of these mixtures. These phases include charge/spin density 
wave order, singlet and triplet pairing, polaron pairing^, and a supersolid phase. Various methods 
to detect our results experimentally are discussed. 



PACS numbers: 03.75.Ss,03.75.Mn,05.30.Fk,05.30. Jp 



I. INTRODUCTION 

Recent advances in controlling ultra cold atoms lead to 
the realization of truly one dimensional systems, and the 
study of many-body effects therein. Important bench- 
marks, such as the Tonks-Girardeau gaa^'^ and the Mott 
transition in one dimension^, have been achieved by trap- 
ping bosonic atoms in tight tubes formed by an optical 
lattice potential. Novel transport properties of one di- 
mensional lattice bosons have been studied using these 
techniques^. More recently, a strongly interacting one 
dimensional Fermi gas was realized using similar trap- 
ping methods^. Interactions between the fermion atoms 
were controlled by tuning a Feshbach resonance in these 
experiments. On the theory side, numerous proposals 
were given for realizing a variety of different phases in 
ultra cold Fermi systems^, Bose- Fermi mixtures-i^iiSiii, 
as well as Bose-Bose mixtureai^. In [l^, commensurate 
mixtures in higher dimensions were studied. 

In this paper we explore the behavior of ultracold 
atomic mixtures, confined to one-dimensional (ID) mo- 
tion in an optical lattice, that exhibit different types of 
commensurability, by which we mean that the atomic 
densities and/or the inverse lattice spacing have an in- 
teger ratio. Commensurable fillings arise naturally in 
many ultracold atom systems, because the external trap 
potential approximately corresponds to a sweep of the 
chemical potential through the phase diagram, and there- 
fore passes through points of commensurability. At these 
points the system can develop an energy gap, which fixes 
the density commensurability over a spatially extended 
volume. This was demonstrated in the celebrated Mott 
insulator experiment by Greiner et al.^^, where Mott 
phases with integer filling occurred in shell-shaped re- 
gions in the atom trap. These gapped phases gave rise to 
the well-known signature in the time-of-fiight images^^, 
and triggered the endeavor of 'engineering' many-body 
states in optical lattices. Further examples include the 
recently created density-imbalanced fermion mixturesJ'^ 
in which the development of a balanced, i.e. commensu- 
rate, mixture at the center of the trap is observed. 

In ID, this phenomenon is of particular importance, 
because it is the only effect that can lead to the open- 
ing of a gap, for a system with short-range interactions. 



In contrast to higher dimensional systems, where, for in- 
stance, pairing can lead to a state with an energy gap, in 
ID only discrete symmetries can be broken, due to the 
importance of fluctuations. Orders that correspond to a 
continuous symmetry can, at most, develop quasi long 
range order (QLRO), which refers to a state in which an 
order parameter 0{x) has a correlation function with al- 
gebraic scaling, {0{x)0{fi)) ~ \x\~^'^~°'\ with a positive 
scaling exponent a. 

Due to its importance in solid state physics, the 
most thoroughly studied commensurate ID system is the 
SU(2) symmetric system of spin-1/2 fcrmions. This sys- 
tem develops a spin gap for attractive interaction and 
remains gapless for repulsive interaction, as can be seen 
from a second order RG calculation. However, the as- 
sumed symmetry between the two internal spin states, 
which is natural in solid state systems, does not generi- 
cally occur in Fermi-Fermi mixtures (FFMs) of ultra-cold 
atoms, where the 'spin' states are in fact different hyper- 
fine states of the atoms. An analysis of the generic system 
is therefore highly called for. Furthermore, we will extend 
this analysis to both Bose-Fermi (BFMs) and Bose-Bose 
mixtures (BBMs), as well as to the dual commensura- 
bility, in which the charge field, and not the spin field, 
exhibits commensurate filling, as will be explained below. 

The main results of this paper are the phase diagrams 
shown in Fig. [THU We find that both attractive and re- 
pulsive interactions can open an energy gap. For FFMs 
the entire phase diagram is gapped, except for the re- 
pulsive SU(2) symmetric regime (cp. Q), for BFMs 
or BBMs the bosonic liquid(s) nced(s) to be close to 
the hardcore limit, otherwise the system remains gap- 
less. Furthermore, we find a rich structure of quasi- 
phases, including charge and spin density wave order 
(CDW, SDW), singlet and triplet pairing (SS, TS), po- 
laron pairing^' and a supersolid phase, which is the 
first example of a supersolid phase in ID. These results 
are derived within a Luttinger liquid (LL) description, 
which treats bosonic and fermionic liquids on equal foot- 
ing. 

This paper is organized as follows: In Section |TT] we 
classify the different types of commensurate mixtures 
that can occur, and in Section IIIII we discuss the effec- 
tive action of the mixtures with the most relevant com- 
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mensurability term. In Section IIVI we discuss the set of 
renormaUzation group equations for such systems, and in 
Section |Vl IVH and IVIIi we apply these results to Fermi- 
Fermi, Bosc-Bosc, and Bose-Fermi mixtures, respectively. 
In Section lVIIIl wc discuss the experimental detectibility, 
and in Section HXl we conclude. 



II. CLASSIFICATION OF COMMENSURATE 
MIXTURES 

We will now classify the types of commensurabil- 
ity that can occur in a system with short-ranged 
density-density interaction. We consider Haldane's 
representatiortlS. of the densities for the two species: 

ni/2 = K/2 + ni/2]^e2™ev2 (1) 

m 

and 1^2 are the densities of the two liquids, 111/2(2:) 
are the low-k parts (i.e. k ^ of the density fluc- 

tuations; the fields Oi/2(a;) are given by Qi/2ix) = 
t:vi/2X + 9i/2{x), with 6'i/2(x) = tt J"^ dylli/2{y). These 
expressions hold for both bosons and fermions. If we 
use this representation in a density-density interaction 
term U12 J dxni{x)n2{x), we generate to lowest order a 
term of the shape U12 J dxlli{x)ll2{x), but in addition 
an infinite number of nonlinear terms, corresponding to 
all harmonics in the representation. However, only the 
terms for which the linear terms (27rmi/2i^i/22;) cancel, 
can drive a phase transition. For a continuous system 
this happens for mivi — m2V2 — 0, whereas for a system 
on a lattice we have the condition miVi — m2V2 = 1^3, 
where mi,m2 and are integer numbers. In general, 
higher integer numbers correspond to terms that are less 
relevant, because the scaling dimension of the non-linear 
term scales quadratically with these integers. We are 
therefore lead to consider small integer ratios between 
the fillings and/or the lattice if present. In [lo!|, we con- 
sidered two cases of commensurabilities: a Mott insula- 
tor transition coupled to an incommensurate liquid, and a 
fermionic liquid at half-filling coupled to an incommensu- 
rate bosonic liquid. In both cases the commensurability 
occurs between one species and the lattice, but does not 
involve the second species. In this paper we consider the 
two most relevant, i.e. lowest order, cases which exhibit 
a commensurability that involves both species. The first 
case is the case of equal filling 1^1 —1^2, the second is the 
case of the total density being unity, i.e. vi + 1/2 = 1, 
where the densities i>i and 1^2 themselves are incommen- 
surate. The first case can drive the system to a spin- 
gapped state, the second to a charge gapped state. We 
will determine in which parameter regime these transi- 
tions occur, and what type of QRLO the system exhibits 
in the vicinity of the transition. These two cases can 
be mapped onto each other via a dual mapping, which 
enables us to study only one case and then infer the re- 
sults for the second by using this mapping. We will write 
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FIG. 1: Phase diagram of a commensurate FFM or a BBM of 
hardcore bosons (with the replacement TSz^SS), in terms of 
the interaction U12 and the parameter z — \vi —V2\/{vi -I-W2). 
For both attractive and repulsive interactions a spin gap 
opens, except for z = and positive interaction. In the at- 
tractive regime, a FFM or a BBM shows either singlet pairing 
or CDW order, or a coexistence of these phases. For repulsive 
interaction these mixtures show SDW ordering, with FFMs 
and BBMs showing subdominant triplet or singlet pairing, 
respectively, for a large range of z. In the gapless regime, a 
FFM shows degenerate SDW and CDW order, and a BBM 
shows SF with subdominant CDW, i.e. supersolid behavior. 
For very large positive values of U12 the system undergoes 
phase separation (PS); for very large negative values it col- 
lapses (CL). 

out our discussion for the case of equal filling and merely 
state the corresponding results for complementary filling. 

III. EFFECTIVE ACTION 

The action of a two-species mixture with equal filling 
in bosonized form is given by: 

S = So^i + So^2 + S12 + Sint- (2) 

The terms S'o.j, with j = 1,2, are given by 

Each of the two types of atoms, regardless of being 
bosonic or fermionic, are characterized by a Luttinger 
parameter K1/2 and a velocity wi/2. Here we integrate 
over r = {vqt,x), where we defined the energy scale 
^'0 — (vi +^2)72. The term 5*12 describes the acoustic 
coupling between the two species, and is bilinear: 

^12 = ^ j d'r 8x9^8^2 + ^ J d\dr0idre2.i4.) 

The second term is created during the RG flow; its pref- 
actor therefore has the initial value Vi2(0) = 0. We define 



3 



So = "So,! + 50,2 + 512, which is the diagonalizable part of 
the action. Sint corresponds to the non-Hnear couphng 
between the two hquids, which we study within an RG 
approach: 



2.912 
(27ra)^ 



J Srcos{29i - 292). 



(5) 



This bosonized description appHes to a BBM, a BFM, 
and a FFM. Depending on which of these mixtures we 
want to describe we either construct bosonic or fermionic 
operators according to Haldane's contructioni^: 



f/b = h + n]i/2 



rn odd/ even 



(6) 



j/q is the zero-mode of the density, ^{x) is the phase field, 
which is the conjugate field of the density fluctuations 
n(a::). The action for a mixture with complementary fill- 
ing, 1/1 + 1/2 = 1, is of the form 5*0 + S'^^f, where the 
interaction S'j'^j is given by: 

2512 



(27ra)2 



(frcos{29i + 2d2 



(7) 



To map the action in Eq. ([2]) onto this system we use the 
mapping: 6*2-^ - 02, </'2^ - 02, and 512^ - 512, which 
evidently maps a mixture with complementary filling and 
attractive (repulsive) interaction and onto a mixture with 
equal filling with repulsive (attractive) interaction. 



IV. RENORMALIZATION GROUP 

To study the action given in Eq. ([2]), we perform an RG 
calculation along the lines of the treatment of the sine- 
Gordon model in [2l|. In our model, a crucial modifica- 
tion arises: the linear combination 9i —02, that appears 
in the non-linear term, is not proportional to an eigen- 
mode of Sq, and therefore the RG flow does not affect 
only one separate sector of the system, as in an SU(2)- 
symmetric system. The RG scheme that we use here 
proceeds as follows: First, we diagonalize So through the 
transformation 



-Bi^l + B202, 
0101+ D202. 



(8) 
(9) 



The coefficients i3i/2 and -D1/2 are given in the Appendix. 
The fields 0i/2 are the eigenmode fields with velocities 
^1/2 (see Appendix). As the next step, we introduce an 
energy cut-off A on the fields 0i/2 according to oj'^/vi/2 + 
Wi/2^^ < A^. We shift this cut-off by an amount dA, and 
correct for this shift up to second order in gi2. At first 
order, only gi2 is affected, its flow equation is given by: 



dgi 



dl 



2-Ki~K2 



2 U12 + V^12t'll^2 
TT Vl + V2 



) 512, (10) 
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FIG. 2: Phase diagram of a BBM with the first species being 
in the hardcore limit, in terms of U12, and the Luttinger pa- 
rameter of the second species {K2), at the fixed velocity ratio 
\vi — V2\/{vi + V2) = 0.5. For large repulsive interaction the 
system undergoes phase separation (PS), for large attractive 
interaction the system collapses (CL). In the regime below the 
thick line the system opens a gap, i.e. if species 2 is close to 
the hardcore limit. However, for larger values of K2, the gap- 
less phase is restored. Close to the transition, the properties 
of the hardcore bosons, are affected by the RG flow, leading 
to supersolid behavior. 



with dl — dA/A. At second order several terms are cre- 
ated that are quadratic in the original fields 0i and 02- 
We undo the diagonalization, Eq. ((8]) and ([9]), and ab- 
sorb these terms into the parameters of the action, which 
concludes the RG step. By iterating this procedure we 
obtain these flow equations at second order in 912: 



dK,/2 


9I2 


dl 




dvi 


9I2 


IT 


lOTT 


dV2 


9I2 


Hi 


— V2 

167r 


dUi2 


9I2, 


dl 




dVi2 


9I2, 


dl 





V2_ 

Vl 



2+ - + - 



"Vi_ 

•~V2 



V2 
V2_ 
Vl 



- V2) 
{l/vi + I/V2) 



(11) 

(12) 

(13) 
(14) 
(15) 



A similar set of equations has been derived in JTj for a 
FFM in non-bosonized form. The difference between our 
result and the result in [2] is the renormalization of the 
velocities, that we find here, which is due to different 
types of expansions: In Q only one- loop contributions 
are taken into account, whereas here we use a cumulant 
expansion in (712, which at second order includes contri- 
butions that are two-loop for the renormalization of the 
velocities. These contributions, which would integrate to 
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zero for equal velocities, as can be seen from Eqns. [12] 
and [131 leads to the discrepancy between the expansion 
in the number of loops and the cumulant expansion, and 
gives a small quantitative correction of the velocities. As 
mentioned before, the advantage of the current approach 
is that the QRLO of the system can be directly deter- 
mined from the resulting renormalized parameters, and 
that the same action can be used to study BBMs and 
BFMs. 

The system of differential equations, Eqns. PH)) to 
([T5)l . can show two types of qualitative behavior: The 
coefficient 1712 of the non-linear term ([5]) can either flow 
to zero, i.e. Sint is irrelevant, or it diverges, leading to the 
formation of an energy gap. In the first case, the system 
flows to a fixed point that is described by a renormal- 
ized diagonalizable action of the type Sq^ from which the 
quasi-phases can be determined. 

When Sint is relevant, we introduce the ficlds^^ 9p/cr — 
^(6*1 ±6*2), which define the charge and the spin sector of 
the system. In this regime, these sectors decouple. Each 
of the two sectors is characterized by a Luttinger param- 
eter and a velocity, Kp^^ and Vp/^, which are related to 
the original parameters in Sq in a straightforward way. 
Using the numerical solution of the flow equations, we 
find that K^r 0, as can be expected for an ordering of 
the nature of a spin gap, leaving Kp the only parameter 
characterizing the QLRO in this phase. 

In order to determine the QLRO in the system we de- 
termine the scaling exponents of various order parame- 
ters. For that purpose, we use the bosonization represen- 
tation of these order parameters, which contain the fields 
O1/2 and 4>i/2, and use the diagonalization, Eqs. ([8]) and 



dH), for the fields 1 
for the fields (j)i/2- 



'1/2 



as well as the dual transformation 



61 = Ci(/)i -I- C202, 

62 = El4'l + £'202- 



(16) 
(17) 



The coefficients C1/2 and Di^2 are given in the Appendix. 
Since the order parameters are now written in terms of 
the eigenfields 61/2 and 1^1/2, the correlation functions 
can be evaluated in a straight forward manner. The scal- 
ing exponents are given by various quadratic expressions 
of the parameters in Eqs. (0, ©, HI]), and HI]). In 
[To| . we give an extensive list of correlation functions, 
which can be transferred to the system considered here, 
with the formal replacement: /3i/2 ^1/2 ; 7i/2 C'i/2j 
Si/2 Di/2, and £1/2 -Ei/2- The order parameter 
with the largest positive scaling exponent shows the dom- 
inant order, whereas other orders with positive exponent 
are subdominant. 



V. FERMI-FERMI MIXTURES 

We will now apply this procedure to the different types 
of mixtures. For a FFM we find that the system always 
develops a gap, with the exception of the repulsive SU(2) 
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FIG. 3: Phase diagram of a BFM with hardcore bosons, in 
terms of the interaction (7i2 and the parameter z = \v\ — 
^^2 1/(^^1 -1-1^2). For both attractive and repulsive interactions 
a spin gap opens, except for 2 = and positive interaction. 
In the attractive regime, a BFM shows either CDW order or 
polaron pairing; for repulsive interaction BFMs show SDW 
ordering. In the gapless regime, a BFM shows CDW order for 
the fermions and SF for the bosons. For very large positive 
values of i7i2 the system undergoes phase separation (PS); for 
very large negative values it collapses (CL). 



symmetric regime (cp. f?']). To determine the QLRO we 
introduce the following operators^^'^*': 



3-cr' 



OsS = y^^S-fli^crSg^a'fL, 

a, a' 

OtS = y^,^fR,<jK,cr'fL,3-. 
a. a' 



(18) 
(19) 
(20) 
(21) 



with (7, cr' = 1, 2, (T = 3 — 2(7, and a — x^y,z. In the 
gapless SU(2) symmetric regime, both CDW and SDW 
show QLRO, with both scaling exponents of the form 
o^SDW/CDW — 1 ^ -^^'p^^i which shows that these orders 
are algebraically degenerate. Within the gapped regime 
the scaling exponents of these operators are given by 
ass,TS, = 2 - and acDW,SDW, ^ 2 - Kp. As dis- 
cussed in [205, the sign of gi2 determines whether CDW 
or SDWz, and SS or TSz appears. In Fig. [T] we show 
the phase diagram based on these results. In addition to 
these phases we indicate the appearance of the Wentzel- 
Bardeen instability, shown as phase separation for repul- 
sive interaction and collapse for attractive interaction. 

We will now use the dual mapping to obtain the phase 
diagram of a FFM with complementary filling from Fig. 
[T] Under this mapping, the attractive and repulsive 
regimes are exchanged with the following replacements: 
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CDW^SDW^, SDW.^CDW, SS,TS,^SDW, and 
SDW^SS. Note that the gapless regime is now on the 
attractive side, with degenerate CDW and SS pairing. 



VI. BOSE-BOSE MIXTURES 

For BBMs we proceed in the same way as for FFMs. 
We introduce the following set of order parameters: 



OcDW 


= blbi+blb2, 


(22) 


Oss 


= bib2, 


(23) 


OsDW^ 


- b\bi-blb2, 


(24) 


OsDW^ 


= b\b2 + blbi, 


(25) 


OsDWy 


= ^l{blb2-b%), 


(26) 



and in addition the superfluid (SF) order parameters bi 
and &2- In Fig- [I]we show the phase diagram of a mixture 
of a BBM of hardcore bosons, which is almost identical 
to the one of a FFM. The phase diagram of the mixture 
with complementary filling, as obtained from the dual 
mapping, is also of the same form as its fermionic equiv- 
alent, with the exception of the gapless regime, in which 
BBMs show supersolid behavior (coexistence of SF and 
CDW order), and with the replacement TSz-^SS. 

In Fig. [21 we show the phase diagram of a mixture 
of hardcore bosons (species 1) and bosons in the inter- 
mediate to hardcore regime (species 2). If species 2 is 
sufficiently far away from the hardcore limit, the sys- 
tem remains gapless. However, in the vicinity of the 
transition the scaling exponents of the liquids are af- 
fected by the RG flow. As indicated, the effective scal- 
ing exponent of the hardcore bosons is renormalized to 
a value that is smaller than 1, and therefore we find 
both SF and CDW order, i.e. supersolid behavior. The 
phase diagram of the dual mixture is of the following 
form: the attractive and the repulsive regime are ex- 
changed, and in the gapped phase we again have the map- 
ping: CDW^SDWz, SDW.^CDW, SS^SDW^^y, 
and SDWx^y^SS . The gapless regime is unaffected. 

The paired SF state discussed in [121 corresponds to 
the SS phase discussed here, whereas the dual SDW^^^j^ 
phase that appears for complementary filling corresponds 
to the super-counter-fluid phase described therein. Note 
that here these orders compete with either CDW or 
SDWz order, and only appear as QLRO, not LRO, as 
in higher dimensions. Both of these insights can only be 
gained by the using the LL description and RG that is 
used in this paper. 
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FIG. 4: Phase diagram of a BFM, in terms of C/12, and the 
Luttinger parameter of the second species {K2), at the fixed 
velocity ratio — V2\/(vi -\- V2) = 0.5. For large repulsive 
interaction the system undergoes phase separation (PS), for 
large attractive interaction the system collapses (CL). In the 
regime below the thick line the system opens a gap, i.e. if 
species 2 is close to the hardcore limit. However, for larger 
values of K2, the gapless phase is restored. Close to the tran- 
sition, the properties of the fermions, are still affected by the 
RG flow, leading to CDW order. 



(see |9lll0l|). and b can develop QLRO in the gapless 
regime. In the gapped regime, the order parameters 



Opp = fpbfLb, 
Opp, = fp.b^fLb\ 



(28) 
(29) 



in addition to Ocdw, show QLRO. {Opp/ppi are special 
cases of the polaron pairing operator (|27p . extensively 
discussed in Q and [lOl-) In Fig- [31 we show the phase 
diagram of a BFM with hardcore bosons, and in Fig. 
[31 we vary the Luttinger parameter of the bosons. In 
both the gapless phase and the gapped phase, we find 
that CDW and /-PP or PP, respectively, are mut ually 
exclusive and cover the entire phase diagram, cp. [9|[lo[. 
The dual mapping again maps attractive and repulsive 
regimes onto each other. Within the gapped phase we 
find the mapping CDW^SDW^, SDW,->CDW, and 
PP-^PP' , the gapless regime is unaffected. 



VIII. EXPERIMENTAL DETECTION 



VII. BOSE-FERMI MIXTURES 

For a BFM we find that the order parameters Ocdw, 
OsDW^, the polaron pairing operator 



Of- 



f-pp 



fpfLe 



(27) 



The phase diagrams that have been derived and shown 
in Figs. [iHH are given in terms of the parameters that ap- 
pear in the effective action. With such a field theoretical 
approach we can find the correct qualitative long-range 
behavior, such as the functional form of the correlation 
functions. However, it is also intrinsic to this approach 
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that the effective parameters appearing can only be qual- 
itatively related to the underlying microscopic parame- 
ters. Based on a phase diagram such as Fig. [21 for in- 
stance, the following features for a mixture of bosonic 
atoms with a short-range interaction can be expected: 
If one species is in the hardcore limit, and the other is 
in between an intermediate interaction regime and the 
hardcore limit, then for attractive interaction between 
them, a gapped state can be created, in which there is a 
competition of SS pairing and CDW order. For repulsive 
interaction, and the second species being very close to the 
hardcore regime, one can also expect a gapped phase, in 
which we find SDW^ order. For the intermediate regime 
we expect a supersolid phase. 

Before we conclude, we discuss how the predictions 
presented in this paper could be measured experimen- 
tally. Since the appearance of a gapped state has already 
been demonstrated for the MI-SF transition in ID 
and since it constitutes a significant qualitative change 
in the system, this would be the first feature predicted 
in this paper to look for. As demonstrated in [2^, RF 
spectroscopy can be used to determine the presence and 
the size of an energy gap. To detect the rich structure 
of QLRO the following approaches can be taken: CDW 
order will create additional peaks in TOF images, corre- 
sponding to a wavevector Q = 2kf. As demonstrated and 
pointed out in (2§|, the noise in TOF images allows to 
identify the different rerimes of both gapped and gapless 
phases. As discussed in |9lflol]. a laser stirring experiment 
could determine the onset of CDW order for fermions, or 
the supersolid regime for bosons. 

IX. CONCLUSION 

In conclusion, we have studied mixtures of ultra-cold 
atoms in ID with commensurate filling. We used a Lut- 
tinger liquid description which enables us to study FFMs, 
BFMs, and BBMs in a single approach. We find that 
FFMs are generically gapped for both attractive and re- 
pulsive interactions, whereas for BFMs and BBMs the 
bosons need to be close to the hardcore limit. We find 
a rich structure of quasi-phases in the vicinity of these 
transitions, in particular a supersolid phase for BBMs, 
that occurs close to the hardcore limit. Experimental 
methods to detect the predictions were also discussed. 

We gratefully acknowledge important discussions with 
D.-W. Wang, A.H. Castro Neto, S. Sachdev, T. Gia- 
marchi, E. Demler, and H.-H. Lin. 

APPENDIX A 

Here we give the coefficients that appear in the trans- 
formations ([5]), (O, (dni), and (fTTl) . that map the original 



fields on the eigenfields at each point in the KG flow. The 
coefficients -B1/2 and -D1/2 are given by: 

Bi = /3iCi+/32Ki, = /3iC2+/32K2, (Al) 

Dl = 51C1+S2K1, D2=SiC2+S2K2. (A2) 

The coefficients (3i/2 and 61/2 are given in [olfioj. where 
the indices 'f and 'b' need to replaced by '1' and '2', 
respectively. The other coefficients are given by: 



Cl = ^/vl/vAC0s9, (2 = \/v2/vASmd, (A3) 
Hi = ~\/v\IVa Sin0, K2 = \/v2lVa COS 9 , (A4) 

where the angle 9 is given by: 

tan 261 = -V12 / V^A^i^A^ - (A5) 
where we used: 



iA = va/V^ + 2Vi2VAf3iSi/Tr, (A6) 

Va = + 2Vi2VAf32S2/7r. (A7) 

The velocities vi/2, corresponding to the eigenmodes 
Si/2, are given by 

(A8) 

where V12 is given by 

Vi2^'2Vi2{PiS2 + P2Si)Itt. (A9) 

VA,a are defined as in (9lfl0l|. The coefficients C1/2 and 
Z)i/2, that appear in the dual transformation, Eqs. (fTBl) 
and pT| , are given by: 

Cl = 71771 -I- 72A1, C2 = 71772 +72 A2, (AlO) 
El = ei?7i -f e2Ai, £^2 = ei^/2 + e2A2. (All) 

7i/2 and £1/2 are given in , with 'f ' and 'b' replaced 
by '1' and '2', and 771/2 and A1/2 given by: 

7?i = \/va/vi cos 6*, 7^2 = \/va/v2 sin 61, (A12) 
Al = -\/va/vi sin 61, A2 = ^Jva/ V2Cos9. (A13) 
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